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A3STRACT 


A  mathematical  model  is  presented  for  the  transient  behavior  of  materials 
which  undergo  thermal  decomposition.  The  general  model  accounts  for  a  rate- 
controlled  decomposition  as  well  as  internal  heat  transfer  by  conduction,  con¬ 
vection,  and  radiation.  The  behavior  of  the  model  under  conditions  of  convec¬ 
tive  heating  is  examined  by  means  of  analog  simulation.  Results  of  this  simula¬ 
tion  for  a  charring  material  and  for  Teflon  are  presented  and  the  results  on 
Teflon  compared  with  experimental  data. 
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I.  INTRODUCTION 


Self- regulating  mass  transfer  systems  for  thermal  protection  have  received  a 
great  deal  of  attention  in  recent  years  primarily  as  a  result  of  problems  as¬ 
sociated  with  the  atmospheric  re-entry  of  long-range  ballistic  and  satellite 
vehicles.  The  status  of  our  understanding  of  the  ablation  phenomenon  has  been 
adequately  covered  by  Adams  ^  and  more  recently  by  Hidalgo,  ^  Steg,  ^  and 
Barriault  and  Yos.^ 

The  major  portion  of  the  theoretical  work  which  has  appeared  dealing  with  the 
ablation  process  has  been  concerned  with  the  behavior  of  glassy  materials  ?“>  ^ 

Predictions  based  on  these  treatments  have  been  adequately  confirmed  by  flight 
test.  No  completely  satisfactory  treatment  of  ablation  in  materials  which  under¬ 
go  thermal  decomposition  has  appeared  in  the  literature.  The  treatment  of 
charring  materials  by  Barriault  and  Yos^  is  of  limited  applicability  and  no  cal¬ 
culations  appear  to  have  been  made  using  their  model.  Treatments  of  non¬ 
charring  materials  such  as  polytetrafluoroethylene  (Teflon)  have  either  con¬ 
sidered  the  material  response  as  incidental®  or  have  treated  the  material  as 
a  purely  conductive  problem.  ^ 

The  purpose  of  this  paper  is  to  present  a  mathematical  model  of  the  ablation 
process  applicable  to  both  charring  and  non-charring  materials.  The  tran¬ 
sient  behavior  of  both  types  of  materials  are  examined  briefly  to  indicate  the 
Utility  of  the  model  in  problems  of  plastic  ablation  under  conditions  of  convec¬ 
tive  heating. 
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II.  MATHEMATICAL  MODEL 


When  subjected  to  environmental  conditions  severe  enough  to  cause  material 
removal,  some  plastics  (such  as  phenolic  nylon)  form  a  more  or  less  adherent 
charred  layer  on  the  healed  surface.  This  char  is  the  end  product  of  the 
thermal  degradation  process  which  proceeds  during  the  heating  process.  Some 
plastics  such  as  Teflon  do  not  appear  to  form  a  char  under  any  practical  set  of 
conditions.  For  the  purposes  of  the  present  paper  these  non-charring  materials 
are  considered  to  be  a  limiting  case  of  the  more  general  behavior  exhibited  by 
the  charring  materials. 

Conceptually  the  original  plastic  is  considered  as  a.  hypothetical  mixture  con¬ 
sisting  of  an  irreducible  "char  matrix"  (subscript  2)  and  a  completely  decompos¬ 
able  "active  material"  (subscript  1).  This  active  material  is  assumed  to  under¬ 
go  thermal  decomposition  yielding  only  gaseous  products.  This  decomposition 
is  assumed  to  follow  a  first  order  rate  law  and  the  rate  constant  is  assumed  to 
vary  with  temperature  in  accordance  with  a  simple  Arrhenius  expression. 

Since  the  material  has  continuously  variable  properties  throughout  it  has  not 
been  necessary  to  define  a  "charring  temperature"  as  in  Barriault  and  Yos^ 
or  even  a  char  layer  as  such.  This  feature  of  the  present  model  has  definite 
advantages,  both  conceptual  and  computational,  over  a  model  requiring  char¬ 
ring  temperatures  as  basic  input. 


Consider  a  control  volume  within  a  bar  of  material  whose  surface  is  receeding 
at  an  unsteady  rate  s  (t  ).  Figure  1  shows  the  geometry  and  coordinate  sys¬ 
tem  employed.  Choosing  unit  cross-sectional  area  for  the  bar,  the  volume 
between  the  planes  x  and  (  x  +  dx  )  is  simply  dx  .  Gaseous  decomposition  pro¬ 
ducts  are  produced  in  unit  volume  at  a  rate  («  pj  ).  This  is  the  rate  at  which 
the  active  component  of  the  material  is  being  converted  to  gaseous  products. 

It  is  assumed  that  the  gases  thus  formed  are  free  to  move  in  the  direction  of 
negative  x  and  that  thermal  equilibrium  exists  between  the  gaseous  and  con¬ 
densed  phases  at  all  points  in  the  material.  Denoting  by  w  the  mass  of  gas 
flowing  across  unit  area  per  unit  time,  continuity  requires  that  equation  (1)  be 
satisfied. 


<9  w 
dx 


P 1 


(1) 


The  total  heat,  dQ  ,  added  to  the  control  volume  in  time  dt  is  the  sum  of  the 
contributions  due  to  thermal  conduction,  dQ.  ,  by  gas  flow,  dQ  .  by  chemi¬ 
cal  reaction,  dQc  ,  and  by  the  absorption  of  radiant  energy,  dQfS 


dQ  a  d  Qjg  +  d  Qg  +  dQc  +  dQ^  . 


(2) 
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Figure  1  GEOMETRY  AND  COORDINATE  SYSTEM 
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Considering  the  one  dimensional  case,  these  conti’ibutions  are  given  by  equa¬ 
tions  (3)  through  (6): 


d  /  dT  \ 

dQ“ =  it  vk  ir)dxdt  ’ 

(3) 

d 

dQg  “  (whg)dxdt 

(4) 

d  Qc  *<  p  j  AHC  dx  dt  , 

(5) 

dF 

8  Qp  “  —  . . dx  dt 

R  dx 

(6) 

In  these  equations,  k  is  the  thermal  conductivity  of  the  bulk  material,  h  is  the 
specific  enthalpy  of  the  gaseous  products,  AHC  is  the  heat  absorbed  per*unit 
mass  of  gases  produced,  and  F  is  the  radiant  flux  at  any  point  in  the  material. 
Assuming  a  constant  specific  heat  for  the  gaseous  decomposition  products  and 
making  use  of  the  thermal  equilibrium  assumption,  dQg  can  be  written  in  the 
form  of  equation  (7); 


dQ 


g 


dxdt 


(7) 


The  heat  added  to  the  control  volume  is  then  given  by  equation  (8), 


do  = 


1* 


0  T 

17 


(?T 


+  CgW_+  MCgT  +  AHC)  -  S'dxdt  ■ 


(8) 


The  heat  equation  is  found  by  equating  dQ  to  the  change  in  enthalpy  of  the  control 
volume  and  neglecting  higher  order  differential  products. 


dt 


<?T 


dx 


dT 


-  *r  'k  ~ ) +  c*w  ~ +  *  +  AHc)  ‘ 


8  dx 


dF 

dx 


(9) 


In  equation  (9)  p  and  h  are  the  density  and  specific  enthalpy  of  the  control  vol¬ 
ume.  Both  the  density  and  enthalpy  are  considered  to  be  additive  functions  of 
the  contributions  due  to  the  "inactive"  and  the  "active"  material  present  in  the 
control  volume.  The  contributions  to  average  density  and  specific  enthalpy  due 
to  the  presence  of  gases  have  been  neglected.  Again,  assuming  constant  speci¬ 
fic  heats  the  density  and  specific  enthalpy  may  be  written  in  the  form  of  equa¬ 
tions  (10)  and  (11); 
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p  -  pj  +  p2 


(10) 


Combining  equations  (10)  and  (11)  yields  equation  (12), 

“  {pCpi  +  P2  (Cp2  -  Cpl)}  T  ,  (12) 


or, 


ph  «  AjpX  +  A2T 


(13) 


A  similar  line  of  reasoning  with  regard  to  the  thermal  conductivity  results  in 
equation  (14), 


k  «  Bjp  +  B2 

where , 


(14) 


B, 


ko-k2 

Po~P2 


(15) 


B1  Po 


(16) 


In  these  equations  the  subscript  o  refers  to  the  original  undecomposed  material. 
The  specific  heat  and  thermal  conductivity  as  defined  by  equations  (13)  and  (14) 
are  not  temperature-dependent  in  the  usual  sense  but  are  instead  functions  of 
the  time-temperature  history  through  their  dependence  upon  the  density  p  , 

The  quantity  is  considered  to  be  a  constant  and  the  quantity  pj  is  assumed 
to  vary  in  accordance  with  a  first  order  decomposition  law. 


Pi  **  -Pi<fi(T) 


(17) 


In  equation  (17)  the  function  $  (  T  )  is  the  first  order  rate  constant  as  a  function 
of  temperature.  This  function  may  be  represented  by  an  Arrhenius  expression 
such  as  equation  (18), 


*(T) 


aeip 


T  +  C 


(18) 
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Combining  equations  (10),  (17),  and  (18)  equation  (19)  results. 


P 


=  (p  -  ?2 )  «  esP 


i 

1  T  +  c  \ 


(19) 


The  radiant  flux,  F  ,  appearing  in  equation  (9)  is  considered  due  to  external 
radiant  heating  and  is  assumed  to  obey  an  exponential  decay  law. 


In  equation  (20)  FQ  is  the  incident  radiant  flux  corrected  for  surface  reflectivity 
and  a  is  an  attenuation  coefficient  which  is  a  function  of  the  temperature,  den¬ 
sity,  and  nature  of  the  material  and  the  wavelength  of  the  radiant  source.  If 
the  incident  radiation  is  not  monochromatic  it  is  necessary  to  use  suitably  wave¬ 
length-averaged  values  for  the  surface  reflectivity  and  attenuation  coefficient. 
The  effect  of  multiple  internal  reflections  is  neglected. 

Using  equation  (20)  the  partial  derivative  appearing  in  equation  (9)  can  be  evalu¬ 
ated, 


BF 

dx 


-»F0  exp 


(21) 


The  set  of  differential  equations  describing  the  system  can  then  be  written  in 
the  following  form. 


(  )  dT  (  ,  B2T  (  Bn)  BT 

)v  +  a2(  -57-  =  jV  +  -£T  +  jCgW  +  B1  irj  *r 


+  | (Cg  -  Aj)  T  +  AHC|  p  +  aF0e*p  |- 


dx|  , 


(22a) 


d  w 

dx 


P 


(22b) 


P 


~(p  -  P2I  a  exP 


f 


T  +  C 


(22c) 
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Equations  (22)  together  with  the  appropriate  initial  and  boundary  conditions  then 
serve  to  determine  the  transient  thermal  behavior  of  an  internally  decomposing 
material. 
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III.  APPLICATION  TO  CONVECTIVE  HEATING 


In  application  of  the  model  to  purely  convective  heating  the  last  term  in  equa¬ 
tion  (22a)  vanishes  since  K0  =  0.  The  initial  and  boundary  conditions  employed 
in  this  analysis  are  as  follows: 


Initial  conditions: 

T(x,  0)  =  constant  =  TQ  r 

(23) 

p  (x,  0)  =  constant  =  p0  • 

(24) 

Boundary  conditions: 

at  *  =  s  (heated  surface), 


Pw' 


|hv  +  t)2  (Ah)| 


+  "x  =  s'fl(Ah> 


(25) 


at  x  =  L  (rear  boundary), 

<">,=L  “  0  » 

(26) 

(21)  . 

(27) 

Equation  (25)  is  a  boundary  condition  which  has  been  discussed  a  number  of 
times  in  the  literature,  ^  and  is  applicable  at  the  stagnation  point  of  a 

hypersonic  vehicle.  In  this  equation  is  the  heat -transfer  rate  calculated  to 
a  cold  isothermal  wall,  Ah  is  the  difference  between  the  enthalpy  of  air  at  stag¬ 
nation  conditions  and  its  enthalpy  at  wall  conditions,  <R  is  the  emissivity  as¬ 
sociated  with  the  material,  is  a  specific  enthalpy  change  associated  with 
surface  recession,  t?.  and  are  transpiration  or  blowing  parameters  associated 
with  the  vapor  species  produced  by  internal  decomposition  and  surface  reces¬ 
sion,  respectively.  These  parameters  can  be  obtained  from  a  knowledge  of  the 
molecular  weight  of  the  vapor  species  and  the  boundary -layer  flow  conditions.  ^ 

In  considering  ablation  process  in  a  finite  bar  of  material  it  is  convenient  to 
introduce  a  transformation  to  a  coordinate  system  based  on  the  receeding  sur¬ 
face.  This  transformation,  first  employed  by  Landau,  is  defined  by  equation  (28), 
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application  of  this  transformation  to  equations  (22)  and  to  the  initial  and  bound¬ 
ary  conditions  yields  the  following  set  of  equations: 


<3T 

IT 


/ 1  y  jBi 

\L  -  s  /  }AlP  +  . 


dr  Uqt-*1vLlAH£(  T,c 


AjP  4  Aj 


(rir)' 


Cgw(L  -  s)  4  (dp/d£)  —  s(L  -  s)  (f  -  1)  (Bjp  4  B2) )  dT_ 


+  A2 


n 


(29) 


5  w 

a? 


c 

-(L-s )(p~p2)ae  T  +  C 


(30) 


dp 

dt 


+ 


dp 

L-s/  <?£ 


( 

-(p-p2)ae  T  +  C 


(31) 


Boundary  conditions : 


at  f  -  0 


{(^ 


<?T^ 


Pw®  |hv  +  i?2  |  +  w£=0,7l(Ah) 


(32) 


at  £  =  1 


Initial  conditions : 
T(£,  0)  =  T0 


0)  =  pQ 


(33) 

(34) 

(35) 

(36) 


In  order  to  completely  determine  the  behavior  of  the  system  an  additional  re¬ 
lation  is  required  between  the  surface  recession  rate,  s  ,  and  one  of  the  other 
parameters  such  as  the  temperature  of  the  receeding  surface,  The  form  of 
this  additional  relationship  depends  on  the  mechanism  of  material  removal  at 
the  surface.  This  mechanism  may  be  due  to  such  effects  as  vaporization, 
combustion,  or  shear  forces.  Relationships  applicable  to  melting  materials 
have  been  considered  by  Fledderman  and  Hurwicz,  ^  and  by  Elotnick  and 
Nordquist,  ®  For  purposes  of  checking  out  the  analog  computations,  simple 
equations  were  written  relating  the  recession  rate  to  the  wall  temperature, 


(37) 
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IV.  COMPUTER  SIMULATION 


Equations  (29)  through  (37)  have  been  programmed  for  analog  simulation  of  the 
thermal  behavior  of  plastic  materials.  As  examples  of  the  two  basic  types  of 
materials  for  which  the  model  was  developed,  calculations  were  made  on  Teflon 
and  on  a  hypothetical  charring  material.  The  calculations  employing  the  char¬ 
ring  material  are  of  interest  only  in  illustrating  the  behavior  of  the  model. 

A.  CHARRING  MATERIAL 

The  properties  employed  in  the  simulation  are  listed  in  Table  1,>  These  pro¬ 
perties  were  used  in  connection  with  a  simulated  satellite  re-entry.  Figure  2 
indicates  the  stagnation  enthalpy  and  cold  wall  heating  employed  for  this  run. 
Figure  3  indicates  the  two  mass  transfer  rates  as  a  function  of  time.  Tempera¬ 
ture  profiles  measured  from  the  original  surface  are  plotted  in  figure  4.  The 
constants  selected  for  this  run  allowed  a  "char"  buildup  of  about  0.  150  inch 
during  the  course  of  the  trajectory.  Several  other  simulations  were  made  using 
constant  heat  inputs.  The  principal  result  of  this  series  of  runs  was  a  demon¬ 
stration  of  the  extreme  transient  nature  of  the  ablation  process  in  decomposing 
materials . 

B.  TEFLON 

A  somewhat  more  detailed  analysis  of  Teflon  ablation  was  undertaken  with  the 
aim  of  verification  of  the  model  by  comparison  with  the  rather  extensive  ex¬ 
perimental  data.  10*  Table  1  lists  the  Teflon  properties  selected  for  use  in 
the  calculations.  The  transpiration  parameter  7j  was  obtained  by  assuming  the 
decomposition  product  to  be  C2F4  and  employing  an  experimental  correlation 
based  on  molecular  weight  for  mass  injection  into  a  laminar  air  boundary 
layer.  1®  The  relation  employed  is  given  by  equation  (38)  where  fi  is  the  mole¬ 
cular  weight  of  the  transpiring  gas, 


The  constants  for  use  in  equation  (3  7)  were  obtained  by  an  analysis  analogous 
to  that  of  Steg.  ^ 

The  simulations  were  made  to  correspond  as  closely  as  possible  to  the  condi¬ 
tions  in  the  experimental  work  of  John  and  Recesso.  11  Figure  5  shows  the 
amount  ablated  as  a  function  of  time  for  several  of  these  runs.  Figure  6  shows 
the  calculated  temperature  distribution  for  one  case  as  compared  to  the  steady- 
state  solution  assuming  constant  properties  and  no  internal  decomposition. 
Surface  temperatures  and  surface  temperature  gradients  computed  during  the 
analog  simulation  agree  with  the  experimental  data  of  Hanst.  12 
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TABLE  I 


MATERIAL  PROPERTIES 


Charring  Material 

Teflon 

P0  (lb/ft3) 

64.4 

137.3 

P2  (lb/ft3) 

32.  2 

0 

Cpl  (Btu/lb- *F) 

0.  34 

0.  30 

Cp2  (Btu/lb-  0  F) 

0.  34 

0 

kQ  (Btu/hr -ft  -  '  F) 

0.  10 

0.  14 

k2  (Btu/hr -ft- °  F) 

1. 00 

0 

C  (Btu/lb  -  “  F) 

0.  20 

0.  20 

a  (sec-3) 

5.  94  x  1012 

4. 70  x  1018 

'  (»R) 

50, 000 

72, 900 

AH.  (Btu/lb) 

500 

680 

fR 

0.  7 

0 

hv  (Btu/lb) 

200 

680 

’Jl 

0.5 

0.43 

^2 

0.5 

0.43 

fi  (in/sec -(°  R)n  i 

1.  96  x  108 

12.  9  x  104 

a 

0 

0.5 

y  CR) 

50, 000 

36,450 
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In  order  to  compare  the  analog  runs  with  the  results  of  experimental  ablation 
measurements,  the  results  of  the  simulations  were  reduced  in  the  same  way 
as  the  experimental  data.  The  results  are  presented  in  figure  7.  The  quantity 
plotted  is  the  ratio  of  the  cold  wall  heat  transfer  rate  to  the  rate  of  mass  loss 
by  the  sample.  The  steady-state  values  predicted  by  the  model  at  two  different 
heat-transfer  rates  are  plotted  for  comparison.  As  can  be  seen  from  figure  7, 
the  transient  nature  of  the  process  as  predicted  by  the  model  shows  excellent 
agreement  with  experiment  when  the  experimental  data  are  interpreted  as 
transient  rather  than  representing  steady-state  behavior. 


-13 


■>) 


0  40  80  120  160  200  240  280 

TIMF  (eornnrtel 

Figure  2  SATELLITE  RE-ENTRY  CONDITIONS 


-14- 


MASS  LOSS  RATE  *  I03 Ob/ ft*- sec) 


80  120  160  200  240  280 

TIME  (seconds) 


Figure  3  MASS  TRANSFER  RATES  FOR  SATELLITE  RE-ENTRY 
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Figure  4  TEMPERATURE  PROFILES  FOR  SATELLITE  RE-ENTRY 
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Figure  5  SURFACE  RECESSION  IN  TEFLON 
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Figure  6  TEMPERATURE  PROFILES  IN  TEFLON 
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Figure  7  TEFLON  HEAT  OF  ABLATION 
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V.  CONCLUSIONS 


The  mathematical  model  presented  for  the  transient  thermal  behavior  of  ma¬ 
terials  which  undergo  thermal  decomposition  is  consistent  with  the  experimental 
data  on  these  types  of  materials.  Preliminary  experience  with  the  model  indi¬ 
cates  the  necessity  for  considering  ablation  in  plastic  materials  as  a  completely 
transient  phenomenon  and  for  intrepretation  of  the  experimental  results  in  terms 
of  this  transient  behavior. 
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